A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
08(08.09.2 SS 433; 09.10.1) 



ASTRONOMY 

AND 
ASTROPHYSICS 



SS433: on the uniqueness of cool relativistic jets 



A. A. Panferov 

Nicolaus Copernicus Astronomical Center, Bartycka 18, Warsaw, PL-00716, Poland; panferov@camk.edu.pl 
Special Astrophysical Observatory, Nizhnij Arkhyz, 357147, Russia; panf@sao.ru 

Received /Accepted 



a ■ 



> ■ 

oo : 

in . 
in ■ 

o : 
o\ . 

Or 
i 1 

o ■ 
a: 



Abstract. The relativistic jets of SS433 are outstanding 
for their optical thermal radiation. The radiation is pro- 
duced by small clouds (10 8 cm) whose lifetime is about 
10 3 times larger than the gas-dynamical crushing time. 
We show that the clouds reside in thermal and dynami- 
cal balance as long as they collisionally interact with the 
wind of the supercritical accretion disk. The interaction is 
caused by the precessional movement of the jets and takes 
place only in the sweep-out zone. Beyond the sweep-out 
zone the interaction ceases and optical jets just terminate. 
The cloud magnetic field amplified in course of movement 
through a medium could play a role in containing a cloud. 
Thus, the clue to the uniqueness of the optical jets of 
SS433 is thought to be their precessional movement, which 
provides an opportunity for collisional interaction of the 
clouds with the wind. 

Key words: Stars: individual: SS433 - ISM: jets and 
outflows 



1. Introduction 

Since the discovery of jets in SS433 in 1978, there has 
been discussion about their uniqueness. By now the zoo 
of various jets has enlarged. There have already been de- 
tected hundreds relativistic jets of AGNs, low velocity jets 
of young stars and a dozen jets of compact stars. Jets of 
SS433 stand out among the others for they harbour cool 
clouds within the relativistic flow. It is not clear, however, 
how the clouds could survive in the relativistic jets and 
be observed. It is also worth mentioning, that the thermal 
fraction is likely to be dominant in the jets of SS433. On 
the contrary, the lack of observational evidence of thermal 
gas in AGN jets strongly restricts its volume-filling factor 
(< 10 -8 ) and implies that it cannot be important in the 
energy budget of these jets ( Celotti et al. 1998 ). These pe- 
culiarities make the jets of SS433 unique. 

Jets of SS433 are highly collimated, 0j < 1°.4, mildly 
relativistic, v- 3 = 7.8T0 9 cms -1 , and powerful, with kinetic 
luminosity ps 10 39 ergs -1 . They appear from funnels 
of the supercritical accretion disk and follow its precession 



rotation with a half opening angle of about 20° and a 
period of P pr = 162 d .5 (see Fig. [j]). At the distance of a 
few 10 12 cm the jets rapidly cool and radiate in the X-ray 
domain. The X-ray emission does not show any structure 
of the jets. Therefore, the X-ray jets are thought to be 
continuous ( [Stewart et al. 1987 ). By contrast, the optical 
jets are observed as fragmented (Vermeulen et al. 1993b). 
Their clumpness may be roughly deduced from energetic 
considerations. The gas must be dense enough to radiate 
the bulk of the emission: 



(1) 



and mass loading of the jet by dense gas must approxi- 
mately satisfy: 



(2) 



s 5 • 10 cm 3 and a volume-filling 
for appropriate parameters: luminos- 



where V ~ 6* 2 i? 3 /4 is the volume of the jet, R is its length, 
n e is electron density, ch q is the H Q emission rate. Thus, 
one finds that n c ? 
factor of / « 10 -5 
ity in H Q line Lu a = 10 36 ergs -1 , total jet kinetic lu- 
minosity Lk = 10 39 ergs -1 , £h q = 10 -24 ergcm 3 s -1 , 
R = R c « 8.4 • 10 14 cm — length of e-folding H Q line 
brightness. It seems that the jets become lumpy in the 
gap between the X-ray part and the optical jet beginning 
at about R in = 2.3 • 10 14 cm ([Borisov k Fabrika 1987|). 



The gas could condense there due to a thermal instability 
as Brinkmann et al. ( Brinkmann et al. 1988| ) have shown 
in a simulation of the thermal evolution of the X-ray jets. 
The instability takes place only in the interior jet, until a 
distance of 10 13 cm. The parameters of the clouds found in 
their work are close to those obtained by Panferov k Fab- 
rika (Panferov & Fabrika 1997) from the Balmer decre- 



ments of the jets: hydrogen density n > 10 cm -3 and 
size of the clouds I < 10 s cm. It should be noted that these 
parameters as well as the filling factor / ks 4 • 10 -6 found 
there are more appropriate for the distance of maximum 
optical brightness of the jets R m ~ 4 ■ 10 14 cm. These pa- 
rameters will be representative ones in what follows. The 
clou ds are grouped in cluster s with sizes of about 10 12 
cm (Panferov k Fabrika 1997). Still bigger structures are 
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quite possible in the jets as may be deduced from the 
structure of the jet emission line. 



Brown et al. (Brown et al. 1995) consider cloudy struc- 
ture of the jets as a result of a temporal switch on and 
off of conditions for thermal instability. They found that 
the parameters of SS 433 are just suitable for the insta- 
bility to occur. However, the instability can work with 
a wide range of the parameters and in SS 433 jets the 
conditions for cloud generation exist rather continuously 
in ti me. Indeed, the data on optical emission of the jet 
(e.g. Vcrmculcn et al. 1993b ) demonstrate that the gen- 
eration of the clouds in the jets is continuous and the 
disappearance of the optical jets happens very rarely. The 
intermittent structure of the jets is natural if the radi- 
ating clouds form in such in situ variable process as the 
thermal instability. The optical jets terminate at a dis- 
tance of about Rj = 3 • 10 15 cm ( [Borisov fc Fabrika 19871) 
and then are observed to be quasi-continuous in ra- 
dio ( Vcrmculcn et al. 1993a ). It is observed as separated 
emission peaks grow and decay at a unchanged wave- 
l ength reflecting an evolution of "bullet" along the jets 
( [Borisov fc Fabrika 1987| ; |Vermeulen et al. 1993b| ). Thus, 
the clouds form in the interior jets and further exist a 
long time giving a signature of the optical jets. 

The actual problem is how do the clouds evolve and 
survive in the relativistic jets? As long as the clouds 
are involved in a motion through a medium, that is, the 
dense wind of the accretion disk, the applied ram pressure 
leads to gas-dynamical instabilities that try to disrupt the 
clouds on a crushing time scale (e.g., Jones et al. 1990| ): 



(3) 



where p is a cloud density, p w is a wind density. For the 
clouds of the jets of SS 433 this time is on the order of 10 2 
s. Besides this, the clouds undergo an evaporation in the 
hotter intercloud medium. The evaporation characteristic 
time may be rather short, too (see Eq. (|l6|) ) . Namely the 
capability of the clouds to survive permits us opportunity 
to observe the jets of SS 433 in the optical. In this paper 
we show that the collisional interaction that seems to be 
the destructive one is a most necessary condition of the 
clouds existence. 

In Sect. 2 we show that the collisional interaction may 
be comparable or even more essential than radiation for 
the thermal state of the clouds. Alth ough this has bee n 
concluded by another authors (e.g., Brown et al. 1995|) , 
we use here the last data on the jets parameters 



( Panferov fc Fabrika 1997 ) and on parameters of UV ra- 
diation of the system — a possible source of heat of the 
clouds ( Dolan et al. 1997 ). In Sect. 3 we discuss evolution 
and confinement of the clouds. 



2. Thermal balance of clouds 

There are two main ways to heat the clouds: by in- 
teraction of jets with the gaseous wind of the super- 
critical accretion disk and by radiation of the system 
( Davidson fc McCray 1980 ). Here we consider competi- 
tion of these two mechanisms. As long as the cloud sound- 
crossing time is shorter than the dynamical time-scale 
<d = r/«j °f the jet, the clouds are in pressure equilib- 
rium with surroundings. The influence of ram pressure of 
the wind on the equilibrium is most essential (see Sect. 3). 
Density in the wind scales as oc r~ 2 , where r is the ra- 
dial distance from the source of the jets. Under plausible 
conditions of approximate constancy of temperature and 
mass of a cloud this implies the following dynamical de- 
pendences of cloud parameters: 



n = 1.6 ■ 10 14 rr, 2 cm" 3 , 



I 



a m7 2/3 

4-10 r 14 cm, 



/ = 1.6 • 10' 



(4) 



where ri4 = r/(10 14 cm) and we have used the values 
of the cloud parameters at i? m , given in Sect.l. In what 
follows, we assume these to be valid everywhere along the 
optical jets. 

Let us compare the rates of cloud heating by the im- 
pinging gas and by radiation. The jets sweep out the 
gas of the wind due to the change of the jet direc- 
tion in the course of the precession and nodding mo- 
t ions and a small amplitu de jitters of unclear nature 
( Margon fc Anderson 1989 ). The energy deposition into 



a cloud by crossing protons with energy e p = m p v? /2 de- 
pends on the collisional thickness, measured by the ratio 
cr c = N/N s , where N = nl = 6.4 ■ 10 : 



2i r -4/3 

14 



cm is the 



column density of the cloud, and N s is the column den- 
sity required to stop a proton. The rate of energy loss of 
a high energy proton in an io nized gas is approximately 
( Ginzburg fc Syrovatskii 1964 ): 



r P = 



2.2 • 10~ 8 .f^n / m„kT 
— 1 + 1.21 — ^ 



-3/2 



erg 



s-\ (5) 



for temperature kT <C m c c 2 , kT -C e p . The factor 
K c = 30A/n c n ~ 1, where A is the Coulomb logarithm. 
In an ionized plasma almost all of the lost energy goes into 
heating the electron gas. The stopping column density of 
a proton with energy e p = 32 MeV is: 



nvj^- w 1.4- 10 



!3 cm 3 . 



(6) 



When N < N s , which is satisfied for the entire jet (see 
Eq. ([§)), the cloud heating rate by the fast protons is: 



tJctpn^Vjl 2 — e p n w Vjnl 3 /N s 



(7) 



This is lowest estimate. The magnetic field may be ampli- 
fied to an equipartition value at the surface of a cloud 
moving through a medium with a weak magnetic field 
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(Jones et al. 1996). For assumed in Eq.(^8|) magnetic field 
the gyro-radius of a 32 MeV proton is 3- Hr ru cm. This is 
smaller than the clouds radius, therefore the trajectories 
of protons may be tangled and the heat input boosted. 

The clouds could also be heated by the colli- 
mated radiation from funnels of the accretion disk 



(Arav & Begelman 1993; Panferov & Fabrika 1993) and 
by the disk UV radiation. The heating is realized by 
(^-electrons — high energy electrons broken away from 
an atom due to ionization. Already for ionization degree 
> 0.5, almost all the energy of (5-electrons is utilized in 
the heating of the electron gas of a cloud. We take into 
account radiation only with frequency beyond the Ly- 
man edge. A cloud intercepts radiation flux with density 
jv = ^{R*/r) 2 B Vl where i?* is radius of the radiation 
source, B v is the Planck function, and relativistic factor 
7~ 3 = (1 — -u 2 /c 2 ) 3 / 2 ~ 1 is omitted. The part of the ra- 
diation which is absorbed is 1 — e~ T , where the optical 
thickness is r = a nn I = a (1 — £) n I, the ionized fraction 
is C — n p/( n p + n n)- Then the radiative heating rate of 
one cloud is: 



H r = i 



I dvj v (l-e- T »), 

J l>-\ 



(8) 



where v\ is frequency of Lyman edge. We adopt here 

cm 2 , for photon with energy 



1.73 • 10 



-22-8/3 
fc keV 



e = IkeV • £kc V in the range 13.6 eV - 12.4 keV, for cosmic 
abundances ( Cruddace et al. 1974 ). 

Parameters of the collimated radiation are not well 
known since its orientation to the Earth is unfavourable. 
Here we adopt rather rough estimates for its lumi- 
nosit y Lx ~ 10 39 ergs -1 a nd temperature Tx ~ 1 
keV flArav fc Begelman 1993| ; |Panferov fc Fabrika 1993[ ). 
Then the radius of a source of the collimated radiation, in 
blackbody approach, is Rx ~ 2 ■ 10 7 cm. The supercritical 
accretion disk is a source of near isotropic UV radiation, 
its parameters in blackbody approach are: Tuv ~ 72000 
1.5 ■ 10 12 cm (Dolan et al. 1997). For given pa- 
gives i? r uv /Hf w 3 ■ 10 5 , i.e. 



K, R d 
rameters and r < 1 Eq 



heating by UV radiation of the disk could dominate heat- 
ing by the X-ray collimated radiation. This ratio might be 
bigger because the clouds have rather r > 1 for the UV ra- 
diation, However situation is complicated by the screening 
of the clouds in a jet. 

The cross-over path length of some particle moving in 
the direction of the jet velocity across the jet is determined 
by the curvature of the jet. The motion of the jet caused by 
nodding is more rapid than that caused by precession and 
its rate is <fi w 5.7- 10~ 7 rads -1 (|Borisov fc Fabrika 19871) . 



Using this rate we obtain the cross-over length for a jet 
with fixed pattern If — v-fi-Jcj) — 3.4 • 10 14 cm, with accu- 
racy of factor 2. Column density along this path is: 



rr+h 



N f = 



dr f n 



1.3- 10 2 



(ru + ku)' 



(9) 



Consequently, we have Nf < N s everywhere in the opti- 
cal jets and the jets are transparent for the rapid protons. 
So, the screening effect is not important for the collisional 
heating. Meanwhile the screening factor of clouds in the 
jet Nf/N may be essential for the UV radiation. We now 
determine Ni = l/a u (l — £) as the column density of a 
cloud layer where r = 1 for the UV radiation. Then the 
screening factor for the UV radiation is Nf/N\ ~ 10 6 (1 — £) 
and essentially depends on the ionization degree. There- 
fore we will consider heating H^ v of the cloud by UV 
radiation as an upper limit of radiative heating (||). Re- 
ally, this heating is bigger than mean heating rate of one 
cloud, on average over the whole jet, in Nf/Ni times; the 
opacity of the jet to UV radiation causes a dominance of 
the collimated radiation in the heating of screened clouds. 

The ratio of the heating rate from fast protons (Eq.(^|)) 
to the radiation heating rate (Eq.fla)) is: 



#uv nN„B*L duB™ 



lOr 



-4/3 
14 l 



(10) 



where we have used the reduced form of Eq. (g) for t > 1 , 



substituting 1 instead of 1 - 
of the accretion disk: 



and used density of wind 



Mr, 



4 7r m n r I Vs K 



1.5 • 10 8 



-3 



(11) 



The mass loss rate and wind velocity are taken to be M w = 
10 -4 M Q /y ( |van den Hcuvcl 1981| ) and v w = 2000 km/s 
(from the width of the emission lines) respectively. Thus 
the heating of the jet by fast protons is comparable at least 
with the radiation heating. Their rates become equal at 
a distance 5.6 • 10 14 cm, until which the main bulk of the 
heat input into the jet takes place and equals to 64% of 
the whole input given by Eq.(13|). The collisional heating 
may be rather much bigger than radiative one, because the 
value of i/ r uv is overestimated approximately in 10 6 (1 — C) 
times due to the screening effect. 

This idea concerning the collisional heating of the jet 
clouds is observationally supported: H Q emission of the 
jets is anisotropic, the maximum of its directional pattern 
is directed to side of jet movement (Panferov et al. 1997). 
This implies that some screening effect exists for the im- 
pinging protons. It is possible if a magnetic field tangles 
trajectories of the protons in a cloud and the stopping 
column density N s becomes smaller than that given by 
Eq.(^|). However, the directional pattern can not be ex- 
plained straightforwardly as a result of a head-on collision 
of the clouds with the wind, because the maximum is in- 
clined to the jet axis. The inclination of the head maxi- 
mum is towards the direction of the precessional motion 
and is about 40° to the jet axis. This rather shows that a 
maximum outcome of the emission line radiation is from 
jet sides moving in the wind. 
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As a result of collision with the wind a cloud deceler- 
ates at a rate: 



Aw 

v 



AM 



At 



(J c n w Vj 
n I 



dr- 



2 ■ 10" 



(12) 



Here we used i? m as the lower limit of the integra- 
tion because evolution of the jet at distances < R m is 
not well established. The deceleration ([l2]) agrees well 
enough with the obse rved value of the deceleration < 10 -2 
( Kopylov ct al. 1987 ). Jets sweep out wind gas with a pat- 
tern speed of <fi t. Then the rate of the heating of the whole 
of optical jet provided by fast protons is: 



H = 



/ dr-^ p n w Or)(0jr) « 2.2 • 10 37 ergs" 1 . 



(13) 



This integral is calculated numerically and dependence of 
Nf on a geometry of the processing conical jet is taken into 
account. The calculated value of the whole heat H agrees 
with the observable radiation output of the jets in H a line 
Ln a ~ 10 36 ergs" 1 flPanferov et al. 1997|): a part < 0.1 of 



absorbed energy is emitted in H Q . 

The small clouds in question would be overheated 
and their H Q radiation would have ceased unless the 
clouds are dense enough to establish a balance between 
radiative losses and the heating by the fast protons. In 
the temperature region near 10 4 K, which is the thresh- 
old ionization temperature of hydrogen, the efficiency of 
the radiative cooling A(T) strongly depends on temper- 
ature and changes by 3 orders of magnitude: A = 10 
- 10 _22 ergcm 3 s" 1 (Kaplan & Pikelner 1979). From this 
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and Eq. (|4|) one finds limits of possible values of the radia- 
tive energy loss of the clouds, corresponding to the end 
and the beginning of the jet: 



3 • 10" 3 erg cm~ 3 s~ 4 < An 2 < 9 • 10 4 erg cm" 



(14) 



Since limits on the heating rate H c from Eq. (0) correspond 
to: 

9 • 10" 2 erg cm" 3 s" 1 < H c /l 3 < 3 • 10 3 erg cm" 3 s"\ (15) 

for the jet end and beginning respectively, the heating rate 
never can exceed the cooling capability, given by inequal- 
ity (0), and the thermal balance of the clouds can be 
maintained near a temperature of 10 4 K throughout the 
whole of the optical jets. 

3. Evolution and confinement of clouds 

The lifetime of the clouds in the jets, about 4 days, is 
much longer than the cloud sound-crossing time. There- 
fore, the clouds must be in pressure equilibrium with the 
ambient medium and should not be exposed to destruc- 
tive processes. The pressure of the ionizing radiation p r = 
H^ v /l 2 c w 9 r^2 dyn is much smaller than the gas pres- 
sure in a cloud p = InkT w 1 3 r^l dyn, where we use a 
temperature of T = 20000 K flPanferov fc Fabrika 19971 ) . 
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Fig. 1. Sketch of precessing jet of SS433. Horizontal axis 
lies in the accretion disk plane. Both axes are in units 
of 10 14 cm. Scales are exact. The wind of the disk fills 
all space except the jet precession surface after distance 
i? S w ~ 3 • 10 15 cm, that is the distance of "crushing wall" . 
And in bounds of the sweep-out zone the wind fills the 
jet precession surface only along a helical pattern moving 
with the precession period. The expansion of the optically 
radiating clouds occurs beyond "crushing wall" between 
two rings. 

Therefore, radiation is not important for dynamical bal- 
ance of a cloud. The pressure equilibrium of clouds with 
surroundings imposes constraints on density = n f ~ 
10 7 cm -3 and on temperature Th = Tn/n^ ~ 10 10 K of 
the surrounding gas at a distance R m . In this hot medium 
the clouds will evaporate with a timescale of: 



pi 



to 6p h c h F(a ) 



« 10 4 



(16) 



where m is the mass of the cloud, Ch is sound veloc- 
ity in the hot medium, F(ao) is some function (given 
by equation (62) in Cowic & McKee 1977), which is ;§> 1 
for our case. The evaporation of the clouds with so short 
timescale contradicts to the lifetime of the clouds. There- 
fore the clouds evaporation must be suppressed. On the 
other hand, the ram pressure of the wind is important in 
pressure balance of the clouds: p w /p — p w v? /2nkT rj 17. 
Ram pressure of the wind gas seems to be strongest and 
must govern the clouds pressure. But again, gas-dynamical 
instabilities are capable of disrupting the clouds on a 
timescale given by Eq.(|). This is the puzzle of SS433 
jets: how can thermal gas survive there? 

The optical jets end at a distance of about 3 • 
10 15 cm. Further, at a distance 3.7 • 10 15 cm, the jets 
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brighten in radio which is a consequence of jet expan- 
sion (Vcrmculcn ct al. 1987). It seems that the expansion 
is connected with the termination of the optical jet. As 
long as the jets precess and move through the dense wind 
of the accretion disk they sweep out the wind gas to the 
surface of the precession cone. After the passage of the jets 
the wind fills the jet channels and during the precessional 
period a region of length R sw = P pl v w « 3 • 10 15 cm can be 
refilled. Consequently, at distances larger than the sweep- 
out distance the lobe along the precession cone surface is 
empty (see Fig. 0) . Thermal diffusion of wind gas is not 
capable of filling the lobe: Ch-P pr <C 9- } r at r > i? sw . The de- 
pendences (^) of the clouds parameters on a distance and 
the clouds heating balance have not peculiarities in point 
r = R sw . But external to the clouds conditions abruptly 
change there. Possibly, this forces the confinement of the 
clouds to switch off and results in the expansion of the jet. 

From laws in Eqs.(Q) we derive the approximate size of 
the cloud I « 4 ■ 10 s cm and the filling factor / w 5 • 10 -7 
at the end of the optical jets. With these parameters the 
clouds start to expand freely and fill the whole volume of 
the jets in time: 



tr 



I 



d .5, 



(17) 



where c s is sound velocity in the cloud. The estimated 
time of expansion corresponds well to the time interval 
required for the jet to move between the end of the op- 
tical part and the radio brightening zone. The last could 
develop due to the increase of the rate of generation of rel- 
ativistic particles, when the clouds expand and turbulence 
is enhanced. 

The coincidence of lengths of the optical jets and the 
sweep-out zone indicates that ram pressure of the wind gas 
causes confinement of the clouds. The ram pressure switch 
off is equivalent to a "hard wall" which crushes the clouds. 
However, the ram pressure by itself is not able to prevent 
clouds from stripping off the gas and thermal evaporation. 
On the contrary, it seems from Eq.(^) that the higher the 
ram pressure the shorter the lifetime of a cloud. In real- 
ity, ionized clouds and their surroundings harbour mag- 
netic fields, which should be important for stability of the 
clouds. The magnetic field inferred from radio observa- 
tions is B « 0.08 G at the distance of radio brightening 



(Vcrmculcn ct al. 1987). Magnetic flux constancy implies 
that B scales as r _1 and we assume magnetic field in the 
jet is: 



B = 3r7} G. 



(18) 



This field could not confine the clouds having internal 

-2 



pressure 



io- 5 



dyn. Jones et al. (Jones et al. 1996) 



showed in their numerical simulations that a "magnetic 
shield" forms around a cloud moving through a medium. 
The field lines of the medium are caught by a cloud, 
stretched and the magnetic field is amplified. The "shield" 
has magnetic pressure comparable to the ram pressure and 



is able to quench gas-dynamical instabilities and prevent 
evaporation of a cloud. Energy of the "magnetic shield" 
is converted from the kinetic energy of the cloud. So, the 
magnetic field needed for the confinement of the clouds 
in the jets of SS433 can be generated by the motion of 
the clouds through the wind of the accretion disk. This 
mechanism of the confinement naturally explains the ter- 
mination of the optical jet beyond the "crushing wall" at 
a distance i? sw , where collisional interaction of the clouds 
with the wind ceases. 

4. Conclusions 

We have considered consequences of collisional interaction 
of clouds in the jets of SS 433 with the wind of the accre- 
tion disk. This interaction appears to be a determining 
factor of the clouds state and evolution. This is empha- 
sized by the fact that the clouds exist only in boundaries 
of the sweep-out zone over which the jet sweeps out the 
wind. We propose that in the sweep-out zone the clouds 
are prevented from destruction essentially by a magnetic 
field which is amplified due to a collisional interaction with 
the wind. The rapid expansion of the clouds after their exit 
from the sweep-out zone naturally explains the disappear- 
ance of the jet hydrogen line emission and the brightening 
of the jets at radio wavelengths. 

The collisional interaction of the jets of SS 433 with its 
surroundings is possible for two reasons: the powerful wind 
from the supercritical accretion disk and the precession 
with cone opening being much larger than the jet open- 
ing angle. These factors are necessary for the existence of 
the cold clouds radiating the hydrogen lines. Their com- 
bination may be the clue to the uniqueness of the optical 
relativistic jets of SS433. 
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